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ABSTRACT 

We present deep neutral hydrogen (H i) observations of the starburst/Seyfert galaxy NGC 3079 
and its environment, obtained with the Westerbork Synthesis Radio Telescope. Our observa¬ 
tions reveal previously unknown components, both in Hi emission and in absorption, that 
show that NGC 3079 is going through a hectic phase in its evolution. The Hi disk appears 
much more extended than previously observed and is morphologically and kinematically lop¬ 
sided on all scales with evidence for strong non-circular motions in the central regions. Our 
data reveal prominent gas streams encircling the entire galaxy suggesting strong interaction 
with its neighbours. A 33-kpc long H i bridge is detected between NGC 3079 and MCG 9-17- 
9, likely caused by ram-pressure stripping of MGC 9-17-9 by the halo of hot gas of NGC 3079. 
The cometary HI tail of the companion NGC 3073, earlier discovered by Irwin et al., extends 
about twice as long in our data, while a shorter, second tail is also found. This tail is likely 
caused by ram-pressure stripping by the strong, starburst driven wind coming from NGC 3079. 
We also detect, in absorption, a nuclear H i outflow extending to velocities well outside what 
expected for gravitational motion. This is likely an atomic counterpart of the well-studied out¬ 
flow of ionised gas present in this galaxy. This may indicate that also large amounts of cold 
gas are blown out of NGC 3079 by the starburst/AGN. Our estimates of the jet energy and 
kinetic power suggest that both the AGN and the starburst in NGC 3079 are powerful enough 
to drive the atomic outflow. 

Key words: galaxies: halos— galaxies: kinematics— galaxies: individual (NGC 3079, 
NGC 3073, MCG 9-17-9)— galaxies: active and starburst — galaxies: jets and outflows 


1 INTRODUCTION 

Deep observations of 21-cm atomic neutral hydrogen (Hi) are 
known to provide a powerful tracer of the structure and evolution 
of galaxies. Given that in general the Hi distribution is often more 
extended than the stellar body of a galaxy, it provides an excel¬ 
lent, and sometimes unique, diagnostic of the role of the environ¬ 
ment and interactions in shaping galaxy properties. It also allows, 
through observations of outflows of cold gas, to study the interplay 
between the ISM and the energetic phenomena associated with ac- 


tive nuclei (e.g.,lGallimore et al.|1999[ Morganti et al.|201: 

and 

star formation (e.g.,|van der Hulst et al.|2001||Boomsma et al. 

2005| 

Walter et al.|2008|lKoribalski et al.|2010|lHeald et al.|2011[ 

Wangl 

et al.|2013>. 


Here we present a detailed study of the Hi in NGC 3079 and 
in its environment. NGC 3079 is well-known for its intense star- 
burst in the central regions and an associated outflow of ionised 
gas, while it also harbours a Seyfert nucleus. Although it has been 
the target of many detailed studies, we have obtained new H i data 


that are much more sensitive than those from previous observations 
with the aim to do a more detailed study of the role of the envi¬ 
ronment in the elevated activity in this galaxy and, vice versa, the 
effect of NGC 3079 on it companion galaxies. 

NGC 3079 is a nearby spiral galaxy which has been studied 
extensively in different wavelengths and resolutions because it is 
known for its several unique properties. One of the most spectacular 
features of NGC 3079 is the ionised gas outflow (super-wind) form¬ 
ing a “super-bubble” at the eastern side of the galactic plane ( |Cecil| 
|et al.|200T] l. The outflow appears to orig inate from the nucl eus and 
reaches velocities up to ~ 1500 kms~^< [Cecil et al.||200l| |. Chan¬ 
dra X-ray observations also show a clear correspondence between 
X-ray emission and the optical Her filaments extending ~ 1.3 kpc 
from the nucleus ( [Cecil et al.|2002| |. Large-scale outflows have also 
been observed in the radio continuum as two, kpc-scale, radio lobes 
extending on both sides of the major axis of the galaxy (jde Bruynj 
[TWTllDuric & Seaquist|19^|Baan & Irwin|1995j l. NGC 3079 is 
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classified as Seyfert 2 ( |Ford|198^ or LINER (Low Ionization Nar¬ 
row Emission Line Region, |Heckman|l 9 80| ). 


The “super-bubble” can be caused by either an AGN (e.g., m 
I win & Seaquist|198^ or a nuclear starburst (e.g., |Sofue et al.|200l| l. 
Indeed, various studies of the nucleus of the galaxy have provided 
evidence for the presence of both a weak AGN and a nuclear star- 
burst. Several X-ray observations provide compelling evidence for 
the existence of an obscured AGN. For instance, the presence of an 
X-ray continuum excess in the 20 - 100 keV band and a 6.4 keV Fe 
line emission ( [Cecil et al.|2002| |, and the detection of both soft 
and hard X-ray emission showing a high h ydrogen column den - 
sity (Ah ~ 10^^ cm"^) towards the nucleus ( [lyomoto et al.||200T i. 
The presence of an AGN is also indirectly supported by Very Long 
Baseline Interferometry (VLBI) observations of [Irwin & Seaquist| 
( [1988[ ) at 5 GHz which revealed aligned multiple-components A, 
B and C (along RA.= 120°) consistent with a pc-scale nuclear ra¬ 
dio jet. This was later supported by several high resolution radio 


continuum observations ( [Baan & Irwin| 19951 [Trotter et al.|1998| 


Sawada-Satoh et al.|2000[ Hagiwara et al.|2004[ Kondratko et al.[ 


2005 [ Middelberg et al.|2007[ l. These high resolution observations 

also revealed additional components (E and F) to the central VLBI 
source with varying spectral indices. Although the additional com¬ 
ponents have similar orientation and spacing as the components de¬ 
tected by l'Irwin & Seaquist| ( [1988[ ), they are not aligned with either 
the kpc-scale radio lobe or the proposed pc-scale jet ( [Kondratko et[ 
[al.[2005[[Middelberg et al.[20F7] l. [Kondratko et alT] ( [2()05[ ) have also 
suggested that the jet may exist with a pc-scale wide angle outfiow, 
which goes to show that the central source is far more complex than 
reported in [Irwin & Seaquist| ( [1988[ ). 

Strong FIR emission ( [Soifer et al.[[1987] l, colour excess 
( [Lawrence et al.[1985] ), and unusually high-density core of molecu¬ 
lar gas relative to other starburst galaxies or AGN (~7200 Mopc"^; 
Planes as et al.[[ 19971 [Sakamoto et al.[[T999[ [Sofue et'ni[2001[ 


Koda et al.[2002[) all p oint towards a nuclear starburst (but see also 
Hawarden et al.|19^ for counter-arguments). 


Previous VLA Hi observations of NGC 3079 by [Irwin ^ 
[Seaquist[ ( [1991[ t showed a smooth disk of H i with however some 
signs of perturbation, e.g. in the southern part. One of the main 
aims of our study is to investigate whether deeper H i observations 
would provide further evidence for such perturbations and whether 
these may be related to the activity observed in NGC 3079. 


[Baan & Irwin[ ( [1995[ ) detected H i and hydroxyl (OH) multi- 
component absorptions towards the centre. The results from these 
observations suggest a combination of rotation and gas outfiow to 
be present in the nuclear region ( [Baan & Irwin| 1995} [Pedlar et al.[ 
\1996\ . The second aim of our work is, exploiting the high spectral 
stability of the Westerbork Synthesis Radio Telescope (WSRT), to 
see whether deeper observations would provide more information 
about this outflow of cold gas. 


NGC 3079 is seen edge-on (/ = 84°; [Heckman|1980'] ) and is 
at a distance of ~ 16 Mpc (i.e. 1 arcsec = 85 pc). NGC 3079 has 
two known nearby companion galaxies, NGC 3073 (dwarf SO) and 
MCG 9-17-9 (small Sb-Sc). The optical appearances of both com¬ 
panions do not show any disturbance. However in Hi, NGC 3073 
appears to be affected by the nuclear activity in NGC 3079 ( [Irwin[ 
[et al.|1987] ). This is observed as an elongated Hi tail coming from 
NGC 3073 pointing away from the nucleus of NGC 3079 in the data 
from [Irwin et al.[ ( [T987[ ). Interestingly, as [Irwin et al.[ ( [1987] ) argue, 
this tail may be caused by ram-pressure stripping of H i from NGC 
3073 by the starburst-driven superwind of NGC 3079. Our deeper 
HI data allow us to study this tail in more detail. 


Table 1. Observational Parameters. 


Pointing centre a (J2000) 
Pointing centre S (J2000) 
Date 

Integration time 

Doppler tracking velocity 
Bandwidth 

Number of channels 

Velocity resolution 

10^ 01”" 57.8" 

+55° 40' 47.0" 

4, 5, 9, 12, 28 Jan 2004 

5 X 12 hr 

1135 kms“^ 

20 MHz 

1024 

16.5 kms“^ 

Line cubes 


Synthesised beam (a x 6) 

rms noise 

19" X 23" 

0.13 mJy beam“^ 

29" X 33" 

0.13 mJy beam“^ 

45" X 48" 

0.17 mJy beam“^ 

78" X 80" 

0.23 mJy beam“^ 

Continuum image 


Synthesised beam (a x 6) 

rms noise 

12" X 14" 

0.08 mJy beam“^ 

Peak flux 

310 mJy beam“^ 


Our observations and the results are presented in Sections [^ 
andj^followed by discussion in Sections 4, 5, and 6. 


2 OBSERVATIONS AND DATA REDUCTION 

NGC 3079 was observed with the Westerbork Synthesis Radio 
Telescope (WSRT). Five 12-hr observations were performed on 4, 
5, 9, 12 and 28 Jan 2004. A bandwidth of 20 MHz with 1024 chan¬ 
nels was used while the observations used Doppler tracking at a 
heliocentric velocity of 1135 kms“^ The parameters of the obser¬ 
vations are summarised in Table [T] 

The data were calibrated and reduced using the miriad pack¬ 
age ( [Sault et aL|[1995[ l. Two standard calibrators, 3C 286, and 
3C 48, were used to calibrate the flux scale and the spectral re¬ 
sponse (bandpass). Separate data sets were generated for line and 
for continuum. This was done by fitting a 4-th order polynomial 
to those channels of each visibility spectrum (’uvlin’) that do not 
contain line emission. 

To improve the dynamic range of the images, self-calibration 
was necessary to correctly determine the antenna gains as a func¬ 
tion of time. This was done using the continuum data where iter¬ 
atively the antenna gains and a model of the continuum sources 
was determined, until convergence was achieved. The calibration 
of the telescope obtained in this way was also copied from the con¬ 
tinuum data to the line data. However, by only following this pro¬ 
cedure, the quality of the images was still limited by the presence 
of strong off-axis sources that could not be calibrated well due to 
direction-dependent imaging errors. These off-axis errors are due 
to small pointing/tracking errors and frequency dependent, quasi- 
periodic variation in the primary beam due to standing waves in the 
dishes of the WSRT. Therefore, we have used the technique com¬ 
monly known as “peeling” to correct for these off-axis errors in the 
image, following relatively standard steps (e.g. [Oosterloo|[2011[ l. 
The procedure involves, after subtraction of the central sources 
from the uv data, self-calibration on the remaining problematic off- 
axis sources (in this case including the well-known double quasar 
and first discovered gravitational lens Q 0957-^561 A, B) to de¬ 
rive extra, direction-dependent corrections. The continuum image 
improved significantly after peeling as most image artefacts were 
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Figure 1. Total H i distribution in and around NGC 3079 at 45" x 48" resolution overlaid onto the optical DSS2 image. Contours are at the levels of 0.05, 0.1, 
0.3, 0.5, 1, 3, 5, 10, 30 and 50 xlO^^ cm“^ . The positions of known and possible companion galaxies are marked, see Table|^for details. 


removed while reaching an rms noise level of 0.08 mJy beam ^ 
which is 32% lower than prior to peeling. 

For the line data, cubes were made using different image 
weightings and taperings. The data were binned a factor 2 in fre¬ 
quency after which additional Hanning smoothing was applied re¬ 
sulting in a final velocity resolution of 16.5 km s~^. Parameters de¬ 
scribing the data are listed in Table 

The cubes were cleaned using the Clark algorithm ( |Clark| 
|1980| ). In order to be able to clean to low flux levels, the clean¬ 
ing was done using masks that were defined by smoothing the data 
to twice the spatial resolution and by identifying regions with line 
emission in these smoothed images. The spatial resolution of the 
final cubes ranges between about 20" and 80" (see Table [^. The 
noise levels, per channel, range between 0.13 and 0.23 mJy beam"^. 
We obtained a 5-cr detection limit for the column density, for a line 
width of 16.5 kms“\ of 3.3 x 10^^ cm“^ for the lowest-resolution 
data cube and 1.2 x 10^^ cm“^ for the highest-resolution data cube. 
The detection limit for small Hi clouds is about 6.5 x 10^ M© (5cr 
for a width of 16.5 km s“^). Total intensity Hi maps were made by 
using masks based on a 2-cr cutoff of data cubes where the resolu¬ 
tion was degraded by a factor 2. 


3 RESULTS 

3.1 The HI morphology and kinematics 

Figure [^illustrates the total H i intensity map showing the distribu¬ 
tion of the gas in and around NGC 3079. The Hi column density 
contours are overlaid onto the Digitized Sky Survey (DSS2) optical 
image. Companions detected in H i are also marked. For a summary 
of the basic H i parameters see Tablej^ Figurej^shows the position- 
velocity (p-v) map taken along the major-axis of NGC 3079. The 
map shows a clear pattern of rotation, with the northern side ap¬ 
proaching and south receding. Very striking is the very asymmetric 
kinematics at all scales. This is further discussed in Sec.|4] 

The improved sensitivity of our observations has allowed to 
detect new, faint and kinematically distinct H i features compared to 
the previous observations. The H i total intensity map of NGC 3079 
shows the H i disk to be far more extended and disturbed than pre¬ 
viously observed ( [Irwin & Seaquist|1991| ). The Hi is seen out to 
a projected radius of ~22 kpc. The integrated Hi flux is ~161 Jy 
kms"^ implying a total Hi mass of ~9.6 x 10^ M©. The improved 
sensitivity of the new data highlights the more extended Hi dis¬ 
tribution resulting in ~25% more Hi detected compared to what 
reported by [Irwin & Seaquist| ( |1991| ). This difference is particularly 
due to the much larger extent of the galaxy west of NGC 3079 in 
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Table 2. Hi parameters of the NGC 3079 Group. 


a(J2000) 

(h m s) 

5(J2000) 

(O , 

Name 

(NED) 

Yys 

(kms“^) 

Mhi 

(10* Mo) 

Label 

10 01 57.8 

-f55 40 47 

NGC 3079 

1116 

96.3 

- 

10 00 48.2 

-f55 36 56 

NGC 3073 

1171 

1.86 

- 

10 01 16.8 

-f55 42 58 

MCG 9-17-19 

1289 

3.42 

- 

10 03 10.0 

-f55 35 56 

1100311.18-^553557.6^ 

1050 

0.09 

A 

10 03 32.5 

-^55 31 19 

1100331.69-^553121.F 

1365 

0.17 

B 

10 01 20.7 

-f55 54 21 

[YGK81] 078 

1164 

0.16 

C 

09 59 42.9 

-f55 53 16 

1095940.91-^555317.8'' 

1246 

0.14 

D 

10 04 43.8 

-f55 19 46 

J100444.00-F551943.2'' 

1148 

0.57 

E 


Note: Objects are labeled according to Fig.[^ 

^ SDSS galaxy with a maximum separation of 0.2' but no redshift information available. 
^ SDSS galaxy with a matching redshift. 



Figure 2. Position-velocity slice taken along the major axis (P.A.=166°). 
Contours are 1, 2, 4, 8, 16, 32, 64, 128 and 256 times 0.225 mJy beam“^ 
and similarly for the negative contour levels (also shown as solid contours). 
The NW side of the galaxy approaching. The dashed lines correspond to the 
range of velocities covered by the H i in emission and are marked for refer¬ 
ence in Fig.[^to be compared with the width of the broad H i absorption. 

the new data. Several new H i extensions and plumes of emission 
are also visible to the south-west and north-east. It is interesting to 
note that in the previous data the galaxy appeared more extended to¬ 
wards the east than to the west, while our deeper observations now 
show that the reverse is true and the H i is much more extended to 
the west, in the direction of MCG 9-17-9. 

Figure |A1| shows selected velocity channel maps (every sec¬ 
ond channel) containing Hi emission detected in NGC 3079. Sev¬ 
eral key properties of the Hi distribution and kinematics can be 
inferred from these images. The emission spans over a wide veloc¬ 
ity range from ~850 to 1410 kms“^ with a systemic velocity of 
~1116 kms“^ The channel maps show that most of the extended 
HI is located in the north-west part of the galaxy. Furthermore, the 
Hi disk appears warped on both north and south sides. This can be 
seen at velocities, e.g., v = 952 and v = 1281 kms"^ At lower 
levels (2cr) there is a sign of a trailing stream of gas visible at ve¬ 
locities near 1050 kms“^ This tail-like feature is likely produced 
by an interaction with the companions seen towards south-east of 
NGC 3079 (see Fig.[^. Other streams and counter-rotating Hi are 
also visible velocities between 1116 and 1165 kms“^ From these 


channel maps it is apparent that NGC 3079 is disturbed and has a 
lopsided nature allover the disk. 

There also seems to be an Hi stream associated with the 
galaxy seen north of NGC 3079 (([YGK81] 078; labeled C in Fig. 
[^; visible for velocities v = 1149 and v = 1182 kms"^). Although 
the outer contours in the channel maps look disturbed, the extended 
gas still shares the disk kinematics (e.g., v = 1182 kms"^). The 
more likely source of the disturbance is an interaction between the 
two galaxies (see Fig.[^. Perhaps, the most interesting feature in 
our data is the new detection of large stream of gas encircling the 
galaxy (Figs. p^and [AT] ). This will be discussed in Sec. |3.4| in rela¬ 
tion to the companion galaxies . 


3.2 The HI absorption 

As expected from previous studies (e.g., |Baan & Irwin|1995'] ), Hi 
absorption is detected against the nuclear continuum emission. The 
position-velocity diagram presented in Fig. shows the absorp¬ 
tion as blue contours. The deep part of the absorption profile ex¬ 
hibits a multiple-component structure with a mixture of broad and 
narrow components, consistent with what found by previous stud- 


ies at higher resolution ( 

Baan & Irwin||1995| |Pedlar et al.||1996| 

|Sawada-Satoh et al.|2000 

(. The strongest absorption extends in ve- 


locity from ~850 kms"^ to ~1440 kms"^ The peak optical depth 
is high, T ~ 0.18, even when assuming that the absorption is ob¬ 
served against the entire continuum source in the inner region (i.e. 
the covering factor / = 1). This corresponds to a column density of 
A^hi ~ 2.0 X 10^^ cm for a spin temperature of Tspin = 100 K. 

|Baan & Irwin| ( |1995] ) have studied this absorption using high- 
resolution observations. Their detected H i absorption covers a sim¬ 
ilar velocity range and shows an even higher optical depth, suggest¬ 
ing that the covering factor for our data is not unity and that the 
absorbing screen is concentrated in the nuclear regions (i.e. not in 
front of the diffuse continuum emission which is included in the 
relatively large central beam of our observations). |Baan & Irwin| 
( |1995| ) identified three distinct absorption components associated 
with the nuclear disk, the inner and the outer disk structures. 

Although our observations cannot provide a comparable spa¬ 
tial resolution, the sensitivity and bandwidth of our data have al¬ 
lowed us to explore the presence of H i (in absorption) not part of a 
regularly rotating structure. In addition to the deep absorption cov¬ 
ering the range of rotation velocities of the H i emission, we detect 
a faint, additional component that covers a velocity range between 
~450 km s"^ and 800 km s"^ well outside the velocity range of the 
regularly rotating disk. At the extreme velocities, it is blueshifted 
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Figure 3. A zoom-in plot of the H i absorption. The vertical line represents 
the systemic velocity. The dashed lines represent the velocity range of the 
Hi emission as indicated in Fig.|^ 


by ~600 km relative to the systemic velocity. To give a better 
view of this faint absorption component, a zoom-in of the H i pro¬ 
file at the location of the peak continuum (which has a fiux density 
of 310 mJy beam“^) is shown in Fig.|^ Interestingly, evidence for 
a blueshifted Hi component was also reported by |Gallimore et al.| 
( |1994| ) based on their VLA observations. The kinematics of the gas 
producing this shallow component is not consistent with gas in a 
simple rotating disk, rather it indicates the presence of outflowing 
gas. Since both starburst winds and AGN are known to be able 
to drive outflows, the superwind or the sub-arcsec radio jet might 
power such this outflow in this galaxy. Both mechanisms are dis¬ 
cussed in Sec.|6] 

This broad and shallow component has a maximum depth of 
~0.35 mJy. Broad and shallow absorption components can be the 
result of errors in calibration caused by, for example, temporal in¬ 
stability of the bandpass of the telescope. In order to exclude this 
possibility, the stability of the bandpass was carefully verified by 
investigating whether similar broad features appear in the spectra 
of the calibrators when the bandpass derived from one calibrator 
observed before NGC 3079 was applied to a calibrator observation 
taken after the science observation. In this way we find that the level 
of instrumental variation of the bandpass is a factor ten below what 
would be needed to explain the faint absorption component, which 
implies that the broad absorption is not an instrumental effect. 

We can estimate the parameters of the broad absorption. If 
we assume that this component is observed against the peak of the 
continuum, its optical depth is r ~ 0.0011. Because there is a large 
uncertainty when using higher resolution data due to the variability 
of the core radio flux density reported in some observations (e.g., 
Sawada-Satoh et aL|2000||Kondratko et aL|2005||Middelberg et al.| 
2007> . we only use this value of r to derive a lower limit to the 
column density. The column density is estimated as A^hi = 1.82 x 
10^^ Tspin f rdv cm This gives us as lower limit of A^hi ~ 4x10^^ 
cm (for Tspin = 100 K and covering factor f = 1). It might well 
be that, if the absorbing gas is located close to the centre of NGC 
3079, that the spin temperature of the gas is well above the 100 K 
assumed here and that these column densities have to be increased 
accordingly. 



Figure 4. Continuum image of NGC 3079 (after “peeling”) overlaid onto 
an DSS2 image. Lowest contour level is 0.2 mJy beam“^ with following 
contour levels increasing factors of 2. 


3.3 The radio continuum 

Figure 1^ shows the continuum image of NGC 3079 overlaid onto 
the optical DSS2 red image. The continuum structure ofNGC 3079 
has been well studied at various resolutions and frequen cies (|Duric| 
|& Seaqui^|1988[ |Baan & Irwin||1995[ |Irwin & Saikia||2003| l. At 

high spatial resolution, it shows two radio lobes (extending about 
2 kpc from the nucleus) extending perpendicularly to the plane of 
the galaxy. These relatively large, extended radio lobes, likely con¬ 
nected to the nuclear outflow, are quite unusual for Seyfert galaxies. 
On larger scales, [Irwin & Saikia| ( |2003| ) found that the radio con¬ 
tinuum follows the galaxy disk. This emission is surrounded by a 
radio halo (extending at least 4.8 kpc from the disk) together with 
other extension/loops particularly clear at low frequencies . 

At the resolution of our observation (~1 kpc), we see some 
indication for the nuclear lobes in the form of small extension per¬ 
pendicular to the disk, although these cannot be well separated from 
the extended disk emission. The peak flux density in the centre of 
NGC 3079 is 310 mJy beam"^ (see Fig.|^. 

The disk emission extends to radii similar to what observed 
by [Irwin & Saikia| ( |2003| ). One outcome from our radio continuum 
data is that the halo appears asymmetric, with the emission on the 
W side being more extended than on the E. On the E side, the extent 
is about 5 kpc while on the S-W side it reaches about 7 kpc. 

This extension to the S-W side was also seen in the 1.2 GHz 
image from [Irwin & Saiki'a[ ( [2003[ ), while the remaining emission 
looks smoother (and larger) in our image. This is possibly the result 
of the improved quality of our image thanks to our direction depen¬ 
dent calibration. We do not see evidence for the radio loops seen at 
low frequencies by [Irwin & Saikia[ ( r2003j l. The larger and asymmet¬ 
ric extent of the continuum halo on the W side partly mimics that 
of the the H i emission. No continuum emission is detected, down 
to the level of 0.24 mJy beam"^ (~3cr), from companion galaxies. 
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Figure 5. Position-velocity slice taken along a line through the centres 
of NGC 3079 and MCG 9-17-9 (RA =290°) showing the bridge be¬ 
tween these two galaxies. The figure also shows the Hi extension towards 
J100311.18-t553557.6. Contour levels are -0.3, 0.3, 0.6, 1.2, 2.4, ... mJy 
beam“^ 


3.4 The NGC 3079 Group 

The NGC 3079 group is known to contain a number of smaller 
galaxies. Our data revealed Hi gas associated with seven galax¬ 
ies that are close (spatially and in velocity) to NGC 3079, and 
thus are likely associated with the group (see Fig. where the ob¬ 
jects are labeled). Among these, four have previously known red- 
shifts. All newly identified members of the group are small, gas- 
rich galaxies. The Hi parameters of each galaxy are given in Ta¬ 
ble]^ A velocity dispersion of cr ~ 80 kms"^ is measured for 
the group. Four of the group members (MCG 9-17-9, NGC 3073, 
SDSS J100311.18-f553557.6 (labeled A) and [YGK81] 078 (la¬ 
beled C)) are seen interacting with NGC 3079. Below, we briefiy 
discuss some details of these interactions. 


3.4.1 MCG 9-17-9 

MCG 9-17-9 is a small spiral galaxy located about 6.5 arcmin 
(~33 kpc) north-west of NGC 3079 (see Fig. [^. It is classified 
as Sb-Sc, and seen fairly face on. Very little is known about this 
galaxy other than from observation of the NGC 3079 environment 
(e.g., |Irwin et al.||T987| |Pietsch et al.||1998| ). Bright 
sion (Lx ~ 1.59 X 10^^ erg s"^) has been detected by 
( |1998| ) which was considered rather high for such a small spiral 
galaxy. The H i appeared undisturbed in the interaction of the group 
in |Irwin et al.| ( |198'7] ). The new Hi data, however, reveal that an in¬ 
teraction between MCG 9-17-9 and NGC 3079 is ongoing. This is 
evident by the Hi ‘bridge’ detected between the two galaxies that 
overlaps both spatially (in projection) and in velocity space, sug¬ 
gesting a physical link between the two galaxies. This can be seen 
clearly in Fig. where the p-v slice along the bridge is shown. The 
south end of the bridge joins the north-west part of the H i disk of 
NGC 3079 which also seems disturbed in that region. The other 
extreme of the bridge connects to MCG 9-17-9. Interestingly, no 
perturbations are seen even in very deep optical images (P.-A. Due, 
priv. comm.), suggesting the bridge is not tidal in nature. We dis¬ 
cuss this in more detail below. 

The HI bridge appears to extend to the other 
side of NGC 3079, towards the other small galaxies 
SDSS J100311.18-f553557.6, and SDSS J100331.69-f553121.1 


X-ray emis- 
Pietsch et al. 



Figure 6. Channels maps from the highest resolution cube (19" x 23") 
showing the two Hi tails of NGC 3073. The lowest contour level is 0.25 
mJy beam“^ which is also the increment of the contour levels. The dashed 
contour corresponds to -0.25 mJy beam“^. 


(marked A and B respectively in Fig. [^. While it is pos¬ 
sible that this means that the bridge extends this far, it is 
also possible that this Hi forms another tail, emanating from 
SDSS J100311.18-f553557.6 because this galaxy seems to have an 
Hi tail pointing towards NGC 3079. Overall, the large streams of 
gas encircling NGC 3079 resemble some of the Hi streams found 
in group galaxies, such as Magellanic Stream (e.g.,|Wannier et al.| 
[T9721 ). 

3.4.2 NGC 3073 

NGC 3073 is a dwarf SABO galaxy located about 10 arcmin (~ 50 
kpc) to the west of NGC 3079. The optical spectrum of NGC 3073 
shows a young stellar population and high levels of star formation 
for its kind et al.||200^ [McDermid et al.||2015| ). [Irwin et ^ 
( |1987j ) revealed an elongated H i tail in this galaxy pointing away 
from NGC 3079, but aligned with the nucleus of NGC 3079. The 
total HI map of this galaxy from our data is shown in Fig.[^ The H i 
tail is clearly visible and turns out to be even longer than injlrwinj 
leF^ l fT^ . However, the tail also seems to consist of two streams. 
This can be seen in the channel maps shown in Fig.[^(from 1132 to 
1198 kms"^) and in Fig.[^ The significant curvature of the tail is 
also visible in Fig.|7] The tail extends to a projected length of about 
5 arcmin (25.5 kpc) and it bends at the end giving a “cometary” 
appearance to the galaxy. 

Although we discuss the Hi tail of NGC 3073 in more detail 
below (Sec. |5.2| ), the presence of the secondary branch may indicate 
that ram pressure stripping is the cause of the tail. The two-tailed 
structure of the Hi tail of NGC 3073 shows a fair degree of sim¬ 
ilarity with that of NGC 4330 ( jChung et al.|2007| ) for which it is 
well established that it is caused by ram pressure stripping by the 
hot medium of the Virgo cluster. 


4 THE STRUCTURE OF THE HI DISK 

Our deep Hi observations of NGC 3079 and its group reveal sev¬ 
eral new features, in and around the galaxy, that are clear signatures 
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Figure 7. Integrated Hi image of NGC 3073 showing the full extent of 
the two tails of NGC 3073. Black contour levels are at 3, 5, 10, 20, 50, 
100 xl0^^cm“^ (45" x 48" resolution), and the grey contour level is 1.5 x 
10^^cm“^ (78" X 80" resolution). 



Radius (arcmin) 

Figure 8. Plot illustrating the lopsidedness ofNGC 3079 (through position- 
velocity slices taken along the major axis). One side of the galaxy has been 
flipped, both spatially and in velocity, on top of the other. It is clear that the 
kinematics of the disk is very asymmetric at all radii. 


of the complex and hectic history of this object. The presence of a 
number of tails and streams, as well as the warped and lopsided 
structure and kinematics of NGC 3079, point to the strong effects 
of interactions with the environment. Although the kinematics of 
the Hi disk of NGC 3079 is dominated by regular rotation, it ap¬ 
pears to have gone through a major “shake up” evidenced by its 
lopsidedness and warped structure. Apart from the morphological 
lopsidedness, in particular visible in the outer regions, Fig.j^shows 
that NGC 3079 also has a strong kinematical lopsidedness in the 
bright disk. This is emphasised in Fig. where the kinematics of 
the two sides of NGC 3079 are compared. It is clear that the kine¬ 
matics of the disk is very asymmetric at all radii. The lopsidedness 
is seen over a large range in radius, but is particularly strong in the 
central regions. This non-standard kinematics, with associated non¬ 
circular motions, may well be related to the elevated activity in the 
central regions (elevated star formation and AGN). 

The distribution of the Hi is asymmetric in the north-south as 
well as in the east-west directions, i.e. more extended towards the 
north and west compared to south and east sides of the disk. This 
gives the galaxy a lopsided appearance. As can be seen in Figs. 
and|^ the morphological lopsidedness of the disk appears to begin 
near the edge of the optical disk where the H i surface density starts 
declining, and slowly increases at larger radii. In particular in the 
direction of MCG 9-17-9 the disk has a large extension. 

Figure shows the position-velocity (p-v) diagram taken 
along the major axis. The disk shows a clear pattern of rotation, 
with the northern side approaching and south receding. Apart from 
this overall pattern of rotation, several streams of gas can be seen 
in this figure. In particular in the northern (approaching) side of 
the galaxy, gas is found over a wide range of velocities outside the 
range of the rotating disk. Similarly, on the receding side gas is 
found at anomalous velocities. This likely corresponds to a stream 
of gas at large radii. 

Lopsidedness in the mass and gas distributions is a well known 
phenomenon and affects a significant fraction of galaxies (e.g., 
[Richter & Sancisi||1994[ [Haynes et al.||1998[ [Sancisi et al.[[2008[ l 

and is generally interpreted as a sign of interaction and/or accre¬ 
tion. Similarly, warping is also common among disk galaxies and is 
known to affect many galaxies where the gas disk extends well be¬ 


yond the stellar distribution fGarci a-Ruiz, Sancisi, & Kuijken|2002[ 
[Jozsa et al.|2009[ ». The warp in the outer disk (north and south) ap¬ 
pears fairly mild but detailed modelling will be needed to determine 
the relative importance of the warp. 

Considering the complicated structure of the gas in and 
around the galaxy, the striking bridge between NGC 3079 and 
MCG 9-17-9, as well as the other gas streams in the direc¬ 
tion of the newly discovered companions ([YGK81] 078 and 
SDSS J100311.18-^553557.6), it is likely that the lopsidedness and 
the warp are induced by interactions with companion galaxies via 
tidal encounters. 


5 THE INTERPLAY BETWEEN NGC 3079 AND ITS 
RICH ENVIRONMENT 

In this section, we discuss in some more detail the Hi tails con¬ 
nected to the companions galaxies MCG 9-17-9 and to NGC 3073. 


5.1 MCG 9-17-9: Tidally or Ram-Pressure Stripped? 


The interaction between MCG 9-17-9 and NGC 3079 is traced by 
an HI bridge as illustrated by Fig. The total H i image of MCG 
9-17-9 shows that the H i bridge forms near this galaxy at a column 
density of about 2 x 10^^ cm“^ and a radius of about 3.3 kpc. How¬ 
ever, Fig. also shows that the bridge connects to the distorted 
outer disk of NGC 3079, suggesting that some part of the gas of 
the entire complex originates from the disk of NGC 3079. The two 
galaxies have a relative velocity of 173 kms"^ with MCG 9-17- 
9 being red-shifted with respect to NGC 3079. A velocity change 
of ~100 kms"^ is seen along the bridge. From the size of the 
bridge and the relative velocities, one finds a kinematical timescale 
(R/V) of ~1.9 X lO^yr, giving an indication of the age of the bridge. 
This timescale is not ve ry different from the a ge of the starburst in 
NGC 3079 estimated by Veilleux et al. < [l994' l (i.e. lO^-lO^yr), sug¬ 
gesting that the interaction between MCG 9-17-9 and NGC 3079 
may have played a role in triggering the starburst. Here, we discuss 
the possible origin of the H i bridge. 
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5.1.1 The case for tidal stripping 

HI bridges or streams between galaxies are often associated with 
tidal encounters. The fact that the outer disk of NGC 3079 seems to 
be affected by the passage of MCG 9-17-9 does indicate that tidal 
effect, at some level, play a role. However, very deep optical images 
(P.-A. Due, priv. comm.) do not show any indication for a tidal dis¬ 
tortion of the stellar distribution of MGC 9-17-9. This makes tidal 
effects unlikely to be the main cause of the H i bridge. Moreover, a 
tidal encounter often leads to two arms, one leading and one trail¬ 
ing, while here we observe only one arm. That tidal effects are not 
the dominant origin of the H i bridge is reinforced from an estimate 
of the tidal radius of MGC 9-17-9. The radius r^ beyond which mat¬ 
ter is tidally stripped due to an interaction between two galaxies can 
be estimated using 


where m is the mass of the satellite, M the mass of the main 
galaxy and R the projected separation between the two ( [Binney ^ 
|Tremaine|1987| t. The un-equality is due to the fact that the projected 
separation gives a lower limit to the true separation. 

The ratio of the masses of the two galaxies we estimate from 
their rotation velocities and their relative sizes. For the rotation ve¬ 
locity at the edge of the regular Hi disk of NGC 3079 (~5 arcmin 
= 23.3 kpc radius) we find 200 km s"^ (see Fig.[^ while for MCG 
9-17-9 we find 40 km s"^ at a radius of -0.7 arcmin = 3.3 kpc. This 
gives a mass ratio of 166. The projected separation between the two 
galaxies is 33 kpc. 

Using these values in Eq. 1 gives a tidal radius of rt > 4.8 kpc, 
larger than the radius of the regular, unaffected disk of MCG 9-17- 
9. However, this estimate is not very much larger than the size of 
the disk of MCG 9-17-9 and it might well be that, e.g. during those 
phases of the encounter between where the satellite may have been 
closer to NGC 3079, tidal effects have played some role. The dis¬ 
tortion of the outer disk of NGC 3079 also points to this. However, 
given that the bridge connects to MCG 9-17-9, the formation of the 
bridge seems to be still ongoing and, given the estimated tidal ra¬ 
dius, it is unlikely that tidal effects are the main reason for the H i 
stream near MCG 9-17-9. 


5.1.2 The case for ram pressure 

Another scenario is that the H i bridge is caused by ram pressure 
stripping of Hi from MCG 9-17-9 by the halo of hot gas NGC 
3079 is known to posses ( [Strickland et al.|2004] ). 

The gas in a galactic disk is stripped if the the ram pressure 
of the hot halo gas is greater than the restoring gravitational force 
per unit area of the galaxy’s disk ( [Gunn & Gott|1972) . The density 
of the hot gas, ph, required for stripping gas at a disk radius R can 
be estimated, assuming a gas disk sitting in a uniform spherical 
stellar/dark matter system, from 

( 2 ) 

where Eg is surface mass density of the gas disk, V^ot the rotational 
velocity, V the velocity of the galaxy through the hot halo, and R 
the radius of the disk under consideration ( [Gunn & Gott|1972T l. 

From the total H i image, we see that the disk of MCG 9-17-9 
appears disturbed for column densities below Eg = 2 x 10^*^ cm"^ at 
a radius of 3.3 kpc. For the rotation velocity we use, as before, 40 
km s"^ while for the velocity V of the galaxy w.r.t. the hot medium 
we assume 200 km s“^ somewhat larger than the observed velocity 


difference between NGC 3079 and MCG 9-17-9. This gives us an 
estimate for the density of the hot gas required for ram pressure 
stripping to occur at the column density level and radius observed 
of 7.9 X 10""^ cm"^. 

This value should be compared with estimates of the density 
of the hot halo of NGC 3079. An extensive X-ray halo, with a lumi¬ 
nosity of Lx - 6 - 8 X10^^ erg s"^ and a temperature of - 4 x 10^ K, 
has been detected around NGC 3079 (ROSAT; [Pietsch et al.|1998| 
Chandra; [Cecil et al.|2002[|Strickland et al.|2004| |. From the param¬ 
eters derived for an exponential distribution of the hot gas (well-fit) 
by [Strickland et al.| ( |2004j ), we get, for the location of MCG 9-17-9, 
an observed density of - 7.3 X 10 cm where fy is a 

volume filling factor. A very similar density was derived by [Pietsch| 
[etal.[ ( [199^ for the hot halo of NGC 3079 assuming a simple spher¬ 
ical distribution. These densities are quite close to that required for 
ram pressure to occurring at the column density and radius we ob¬ 
serve the HI disk of MCG 9-17-9 to be affected. 

From the above analysis of both tidal and ram-pressure strip¬ 
ping, we conclude that the stripping that is currently occurring for 
MCG 9-17-9 is mainly due to the interaction with the halo of hot 
gas known to surround NGC 3079. The observed densities of the 
halo match quite well the estimated required densities for the strip¬ 
ping to occur at the observed column densities and radius. The esti¬ 
mated tidal radius is larger than the size of the disk of MCG 9-17-9. 
However, the difference is not very large and it is conceivable that, 
during an earlier phase of the encounter where possibly MCG 9-17- 
9 was closer to NGC 3079, tidal effects may have been important. 
The distorted outer H i disk of NGC 3079 also suggest that this is 
the case. 

5.2 The two tails of NGC 3073 

The ‘cometary’ appearance of the Hi in NGC 3073 was first ob¬ 
served and discussed by [Irwin et al.[ ( [1987] ). Our deeper data reveal 
that the Hi tail of NGC 3073 is not only more extended, but also 
that it appears to be composed of two tails superimposed (Fig.[^. 

Irwin et al. discuss the interesting possibility that this tail 
of NGC 3073 is due to ram-pressure stripping by the fast star- 
formation driven wind coming from the central regions of NGC 
3079. This possibility was also discussed by [Strickland et al.[ 
( [2004[ ). Before we discuss this option in more detail, we first con¬ 
sider the alternatives of tidal effects and ram-pressure stripping by 
the hot halo of NGC 3079. 

The first option can likely be excluded. In the previous sec¬ 
tion we showed that tidal effects are currently not very strong for 
MCG 9-17-9 and there are a number of reasons why for NGC 3073 
this is even more the case. In the first place because the distance 
between NGC 3079 and NGC 3073 is likely to be larger, given 
the larger projected separation of 50 kpc. Moreover, NGC 3073 is 
more massive than MCG 9-17-9, reducing the tidal impact of NGC 
3079. Using the results from [Cappellari et ^ ( [2013[ ), we estimate 
the rotation velocity of NGC 3073 to be 86 km s"^ compared to 40 
km s"^ for MCG 9-17-9. Taking this together, we find that the tidal 
radius for NGC 3073 is likely to be larger than about 17 kpc, which 
is much larger than the size of the galaxy (~1.2 kpc radius), ruling 
out a tidal origin for the H i tail. 

Similarly, we can also rule out that the general halo of hot gas 
is stripping the Hi from NGC 3073. In an similar way as done for 
MCG 9-17-9, we estimate that the required density of the hot halo 
in order to be able to strip the gas from NGC 3073 is 1.3 x 10“^ 
cm“^ which is much larger than the observed densities. 

However, as [Irwin et al.[ ( [1987] ) remark, it is possible that the 
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ram pressure from the wind coming from NGC 3079 is stripping 
NGC 3073. This is mainly, since ram pressure depends on ve locity 
as v^, because of the high velocity of the wind (1000 kms“^ Cecil 


|et al.|200T] ). For this velocity, the required density to strip the H i at 
the column density observed (2.5 x 10^^ cm“^), is about 7.5 x 10“"^ 
cm“^. The lower limit for the electron density of the super-wind, 
near the starburst in the centre of NGC 3079, is estimated to be 
Hq = 2.7 cm"^ (Strickland et al. 2004). Assuming that the density 
of the wind declines with distance from the starburst according to 
r"^, we obtain a density of about 1 x 10“^cm“^ at the location of 
NGC 3073. Given the uncertainties, this matches the required den¬ 
sity for stripping and we reiterate the suggestion made by |Irwin| 
|et al.|(T987] ) that the cometary H i tail of NGC 3073 is due to ram- 
pressure stripping by the fast, starburst driven wind from the centre 
ofNGC 3079. 


6 THE INTERPLAY BETWEEN THE CENTRAL 
ACTIVITY AND THE Hi DISK 

The finding of the H i outfiow seen in absorption can be added to 
the many interesting features of the nuclear region of NGC 3079. 
The broad and shallow blueshifted Hi absorption component is 
blueshifted up to 600 kms"^ with respect to the systemic velocity. 
Being detected in absorption, this ensures that the gas is in front 
of the radio continuum and the blueshifted velocities indicate that 
the absorbing gas is outflowing. Such a broad, blueshifted compo¬ 
nent has been tentatively seen also in the VLA observations of |Gal-| 
[limore et~^ ( |1994| ). A component with such large velocities must 
be associated to extreme, non-gravitational radial motions. 

Fast Hi outflows have been detected in several other radio 
sources with a range of AGN power. Many cases detected so far 
have been found to be associated to powerful radio jets (e.g., |M or^ 
Iganti et ar]|2005| 2013), although cases associated to low-power 
radio sources (e.g., NGC 1266; |Nyland et al.|2013| Mrk 231; |Mor-] 
|ganti|20lT| ) are also known. In a number of these objects, the most 
likely origin for the outflow is considered to be the kinetic push 
by the radio plasma jet following its interaction with a dense ISM 
(e.g. |Morganti et al.|2013| t. This has also been suggested to occur in 
lower-power radio sources such as NGC 1266 ( |Nyland et al.|2013| l, 
a source with a radio power more similar to that of NGC 3079. 

In the case of NGC 3079, it is hard to identify the origin of 
the outflow. This is due to the difficulty to disentangle the effects 
due to the starburst from those caused by the AGN, as it has been 
discussed in detail for the outflow of ionised gas by |Veilleux et al.| 
^l994\ and |Cecil et al.| ( |200^ . According to these studies, the star- 
burst component is the most likely driver, although an AGN com¬ 
ponent cannot be excluded. Indeed, while the high infrared lumi¬ 
nosity and the morphology of the well-known bubble observed in 
NGC 3079 suggest the stellar wind to be important( |Veilleux et al.| 
|1994| ), it has been suggested that the energy output of the AGN and 
starburst activities are comparable ( [lyomoto et akpOOT) . 

The possibility of the existence of a molecular outflow driven 
by the active nucleus was already suggested by [Hawarden et al.| 
( |1995| ). Molecular outflows have now been observed in a number 
of objects. The multi-phase character of gas outflows, and the co¬ 
presence of molecular and atomic hydrogen (as well as ionized 
gas), has been established in a growing number of cases where 
the outflows are driven by AGN. Thus, the Hi outflow detected 
in our observations may represent the counterpart of the outflow¬ 
ing molecular gas (e.g. as in the case of the radio-loud Seyfert 2 
IC 5063, |Morganti et al.|2015[|Tadhunter et al.|2014| l. 


The limited spatial resolution of our observations cannot add 
information about the precise location of the outflow. However, in 
order to shed light on its origin, we can at least derive the energetics 
of the outflow and compare this with what is associated with the 
starburst and the non-thermal radio emission. This should tell us 
whether any of them can be ruled out because of insufficient energy. 

From the lower and upper limit of the column density as de¬ 
scribed in Sec. |3.2| a range of values for the mass outflow rate can 
be obtained following |Heckman] \2002\ and |Rupke et al.| ( [200^ : 

Vhi 


M~30|-^ 


(S)(T^)(To^)(30o^)'^°y’"‘ 

where the gas is outflowing at a velocity v from radius n over a 
solid angle Q. If we assume the outflow originates from the nuclear 
region, r = 100 pc, Q to be :7r steradians and v = 500 kms"^ the 
above equation yields a lower limit to the mass outflow rate of M 
ranging from 0.05 to 0.5 M© yr"^ where we have allowed to 
range between 100 K and 1000 K. These mass outflow rates are 
low compared t hose derived for ULIRGs (ranging between 10 to 
1000 Mq yr"^ ( Heckman]|2002 Rupke et al. 2002 Rupke et al. 
2005 a|b| and for radio gal axies, ranging between 1 and 50 M© yr"^ 


(e.g. 


Morganti et al.|2005] ). 

Veilleux et al.| ( |1994| derived an upper limit for the ionised gas 


outflow in NGC 3079 of M < 10 x n~^ Mq yr"^ (for 10 < < 125 

cm"^). Thus, the ionised gas outflow rate in NGC 3079 is relatively 
modest compared to ULIRGs, but of similar amplitude as the H i 
outflow rates. 

With the Hi mass outflow rates, we can derive the kinetic 
power of the outflow using the approach of |Holt et al.| ( |2006} , in¬ 
cluding both the radial and turbulent component of the gas motion. 


M 

E ~ 6.34 X 10"'’ I — 


FWHM" 

1.85 


erg s 


(4) 


The estimated values range between a few xlO^^ and 10"^^ erg s"^ 
again reflecting the range of column densities. 

The energy associa ted with the nuclear starburst as derived 
by Veilleux et al. (|1994| ) is estimated to be ~6 x ergs 


Veilleux et al.| { |1994| ) quote that for values of > 40 cm the 


injection rate of the nuclear starburst would have been able to pro¬ 
duce the Hi outflow (i.e. providing a power of > 10"^^ erg s"^). 

Recent results (e.g. Nyland et al. 2013) have shown that fast 
HI (and molecular gas) outflows can also be driven by low-power 
radio sources due to the high efficiency of the coupling between 
the radio plasma with the surrounding ISM/IGM e.g.|Cavagnolo et| 
l^ pOTO] ). Thus, it is worth exploring the energetics connected to 
the radio plasma also for NGC 3079. Indeed, the presence of an 
interaction of the parsec-scale jet with ambient gas in NGC 3079 
has been shown by |Cecil et al.| ( |200l'| ) and |Middelberg et aL] ( |2007| ). 

A number of studies have recently addressed the crucial is¬ 
sue of the conversion between radio luminosities and jet power. 
We use the relations proposed by |Willott et al.| ( p^^ , |Wu| ( |2009| ), 
[Birzan et~^ ( |2008| ) and [Cavagnolo et al.| ( |2010| l. Despite the dif¬ 
ferent approaches used by these authors in deriving these relations, 
we obtain a reasonable agreement from these different methods, in¬ 
dicating a jet power of ~4 x 10"^^ erg s"^ Thus, also the jet power 
appears to be large enough to drive the H i outflow. This is further 
supported by the presence of jet-cloud interactions in the galaxy, as 
suggested by |Cecil et al.| ( |200l] ) and |Middelberg et ar] ( |2007| ). 

However, a caveat to the radio plasma driving the outflow is 
whether the putative VLBI jet is really connected to the radio lobes. 
Interestingly, the position angle of the VLBI scale radio jet is close 
to the plane of the galaxy ( [Irwin & Seaquist|1988| l, and misaligned 
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by about 65° with respect to the large-scale lobes. In this situation, 
one could hypothesise that this is due to a deflection partly caused 
by the interaction with the medium: this would be supported by 
the presence of broad blue-shifted Hi absorption. A similar (albeit 
not so extreme) misalignment was found in at least an other spiral 
galaxy - IC 5063 - where the jet is seen in the plane of the galaxy. 
Interestingly, IC 5063 also hosts massive outflows of ionised. Hi 
and molecular gas (|Morganti et al.||2007| 2015; |Tadhunter et al.| 
|20T4] ). 


Thus, our considerations show that it is plausible for the radio 
jet to have a role in driving the H i outflow. However, the connection 
between the complex radio continuum structure in the inner pc re¬ 


gion ( 

[Trotter et al.[1998|[Sawada-Satoh et al. 

2000 [Hagiwara et al.[ 

|2004| 

, Kondratko et al.[2005| Middelberg et a 

.[2007J and the large- 


scale disk, and also the connection between the pc ( [Kondratko et al.| 
|2005| ) and kpc-scale outflows in this galaxy still remains ambigiu- 
ous. 


7 CONCLUSIONS 

We have presented deep Hi observations of the famous star- 
burst/Seyfert galaxy NGC 3079 and its surroundings. The improved 
sensitivity of the data compared to those in the literature have al¬ 
lowed to uncover new and interesting features which give an im¬ 
proved picture of the hectic evolution of NGC 3079 and its envi¬ 
ronment: 

(i) NGC 3079 has a more extended (out to ~22 kpc radius) and 
disturbed H i disk than previously observed. The H i disk is lopsided 
on all scales, both morphologically and kinematically, and shows a 
mild warp on both the north and south sides. The strong kinematical 
lopsidedness in the central regions, and the non-circular motions 
associated with that, are possibly connected to the high levels of 
star formation in the central regions. 

(ii) NGC 3079 is clearly interacting with its neighbours. Large 
streams of gas and tails are visible in the NGC 3079 group. The 
disk of NGC 3079 is strongly affected by these interactions. 

(iii) In addition to the two previously known companion galax¬ 
ies (NGC 3073 and MCG 9-17-9), Hi gas has been detected from 
five more galaxies (three with previously unknown redshift) which 
are likely associated with the NGC 3079 group. Among these, two 
of them ([YGK81] 078 and SDSS 1100311.18-^553557.6) seem to 
be taking part in the interactions in the group. 

(iv) A long HI bridge is visible between NGC 3079 and MCG 9- 
17-9 spanning over a projected length of ~33 kpc. We suggest that 
ram pressure due to the hot X-ray halo of NGC 3079 is most likely 
responsible for the stripping of gas from MCG 9-17-9, although 
also tidal effects may have been relevant in the past. 

(v) The cometary Hi tail of NGC 3073, earlier discovered by 
[Irwin et ^ ( |1987| ), extends even further, while a second, shorter 
tail is now also revealed. Our result agrees with the conclusion of 
[Irwin etar] ( |1987| ) that ram pressure by the super wind coming from 
NGC 3079 is responsible for this tail. 

(vi) An HI counterpart to the ionised gas outflow of NGC 3079 
has been detected in the form of a blueshifted wing of the H i ab¬ 
sorption profile against the central continuum source, extending to 
velocities well outside the range of rotational velocities. However, 
we cannot distinguish between the starburst and AGN driving this 
outflow, with both having comparable power sufficient to drive the 
outflow. 
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Figure Al. Selected channel maps at 45" x 48" resolution, in every second channel, containing Hi emission. Lowest positive contour level 0.3 mJy beam“^ 
which is also the step in contour level. The negative contour is -0.3 mJy beam“^ The channel velocity (in kms“^) is marked at the top left corner of each 
panel. 
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Figure A2. Fig. |Al| continued. 















